scPipe report for sample LC396

Parameters for each preprocessing step

Parameters for sc_trim_barcode

File paths

• input fastq1: /stornext/General/data/user_managed/grpu_mritchie_1/hickey/SCORE/C075_Grant_Coultas/extdata/NN176/scPipe/LC396/LC396.R2.fastq.gz
• input fastq2: /stornext/General/data/user_managed/grpu_mritchie_1/hickey/SCORE/C075_Grant_Coultas/extdata/NN176/scPipe/LC396/LC396.R1.fastq.gz
• output fastq: /stornext/General/data/user_managed/grpu_mritchie_1/hickey/SCORE/C075_Grant_Coultas/extdata/NN176/scPipe/LC396/LC396.combined.fastq.gz

Read structure

assume read1 contains the transcript

• barcode in read1: NA
• barcode in read2: start at position 6, length 7
• UMI in read2: start at position 0, length 6

Read filter

• remove reads that have N in its barcode or UMI: FALSE
• remove reads with low quality: TRUE
  • minimum read quality: 20
  • maximum number of base below minimum read quality: 2

Parameters for alignment

• input fastq: /stornext/General/data/user_managed/grpu_mritchie_1/hickey/SCORE/C075_Grant_Coultas/extdata/NN176/scPipe/LC396/LC396.combined.fastq.gz
• output bam file: /stornext/General/data/user_managed/grpu_mritchie_1/hickey/SCORE/C075_Grant_Coultas/extdata/NN176/scPipe/LC396/LC396.combined.subread.bam
• genome index: /stornext/General/data/user_managed/grpu_mritchie_1/hickey/SCORE/C075_Grant_Coultas/extdata/GRCm38.p6/GRCm38_with_ERCC

Parameters for sc_exon_mapping

• input bam file: /stornext/General/data/user_managed/grpu_mritchie_1/hickey/SCORE/C075_Grant_Coultas/extdata/NN176/scPipe/LC396/LC396.combined.subread.bam
• output bam file: /stornext/General/data/user_managed/grpu_mritchie_1/hickey/SCORE/C075_Grant_Coultas/extdata/NN176/scPipe/LC396/LC396.combined.exon.bam
• transctiptome annotations: /stornext/General/data/user_managed/grpu_mritchie_1/hickey/SCORE/C075_Grant_Coultas/extdata/GRCm38.p6/gencode.vM18.primary_assembly.annotation.gff3, /stornext/Home/data/allstaff/h/hickey/.local/share/renv/cache/v5/R-3.6/x86_64-pc-linux-gnu/scPipe/1.8.0/96dbe230d12477825fb3768b2fa4c3a5/scPipe/extdata/ERCC92_anno.gff3
• do strand specific mapping: TRUE
• fix chromosome names: FALSE

Parameters for sc_demultiplex

• input bam file: /stornext/General/data/user_managed/grpu_mritchie_1/hickey/SCORE/C075_Grant_Coultas/extdata/NN176/scPipe/LC396/LC396.combined.exon.bam
• output folder: /stornext/General/data/user_managed/grpu_mritchie_1/hickey/SCORE/C075_Grant_Coultas/extdata/NN176/scPipe/LC396
• barcode annotation file: /stornext/General/data/user_managed/grpu_mritchie_1/hickey/SCORE/C075_Grant_Coultas/extdata/NN176/scPipe/LC396/LC396.barcode_annotation.csv
• maximum mismatch allowed in barcode: 1

Parameters for sc_gene_counting

• output folder: /stornext/General/data/user_managed/grpu_mritchie_1/hickey/SCORE/C075_Grant_Coultas/extdata/NN176/scPipe/LC396
• barcode annotation file: /stornext/General/data/user_managed/grpu_mritchie_1/hickey/SCORE/C075_Grant_Coultas/extdata/NN176/scPipe/LC396/LC396.barcode_annotation.csv
• UMI correction: simple correction and merge UMI with distance 1
• gene filtering: FALSE

Data summary

The organism is “mmusculus_gene_ensembl”, and gene id type is “ensembl_gene_id”.

Overall barcode statistics

if (is.null(params$organism) || is.na(params$organism)) {
  sce = create_sce_by_dir(params$outdir)
} else {
  sce = create_sce_by_dir(params$outdir, organism=params$organism, gene_id_type=params$gene_id_type)
}
overall_stat = demultiplex_info(sce)
datatable(overall_stat, width=800)

Plot barcode match statistics in pie chart:

plot_demultiplex(sce)

Read alignment statistics

ggplotly(plot_mapping(sce, dataname=params$samplename, percentage = FALSE))

ggplotly(plot_mapping(sce, dataname=params$samplename, percentage = TRUE))

Summary and distributions of QC metrics

if (any(colSums(counts(sce)) == 0)) {
  zero_cells = sum(colSums(counts(sce)) == 0)
  sce = sce[, colSums(counts(sce)) > 0]
} else {
  zero_cells = 0
}

Datatable of all QC metrics:

sce = calculate_QC_metrics(sce)

Warning: 'spikeNames' is deprecated.
See help("Deprecated")

Warning: 'isSpike' is deprecated.
See help("Deprecated")

Warning: 'isSpike' is deprecated.
See help("Deprecated")

Warning: 'isSpike' is deprecated.
See help("Deprecated")

Cache found
Cache found

if(!all(colSums(as.data.frame(QC_metrics(sce)))>0)){
  QC_metrics(sce) = QC_metrics(sce)[, colSums(as.data.frame(QC_metrics(sce)))>0]
}
datatable(as.data.frame(QC_metrics(sce)), width=800, options=list(scrollX= TRUE))

Summary of all QC metrics:

datatable(do.call(cbind, lapply(QC_metrics(sce), summary)), width=800, options=list(scrollX= TRUE))

Number of reads mapped to exon before UMI deduplication VS number of genes detected:

ggplotly(ggplot(as.data.frame(QC_metrics(sce)), aes(x=mapped_to_exon, y=number_of_genes))+geom_point(alpha=0.8))

knitr::knit_exit()